Renormalizat ion- group description of nonequilibrium critical short-time 
relaxation processes: a three-loop approximation 
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The influence of nonequilibrium initial values of the order parameter on its evolution at a crit- 
ical point is described using a renormalization group approach of the field theory. The dynamic 
critical exponent 8' of the short time evolution of a system with an n-component order parameter 
is calculated within a dynamical dissipative model using the method of e-expansion in a three- 
loop approximation. Numerical values of 9' for three-dimensional systems are determined using the 
Pade-Borel method for the summation of asymptotic series. 
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I. INTRODUCTION 

This study is devoted to the influence of nonequi- 
librium initial states on the evolution of magneti- 
zation m(t) of a ferromagnetic system at a crit- 
ical point. As is known [1], anomalous features 
in the phenomena of critical dynamics are deter- 
mined primarily by the long-range correlation of 
long-lived fluctuations of some thermodynamic vari- 
ables. In this context, subject of fundamental inter- 
est to study is the process of critical relaxation of 
a system from an initial nonequilibrium state - for 
example, that created at temperatures much higher 
than the critical temperature and, hence, character- 
ized by a short correlation length - to a strongly 
correlated state at the critical point. Janssen et al. 
[3] showed that the critical evolution of a system 
from the initial nonequilibrium state with a small 
magnetization m$ — m(0) <C 1 displays a univer- 
sal scaling behavior of m(t) over a short time early 
stage of this process, which is characterized by an 
anomalous increase in magnetization with time ac- 
cording to a power law. The exponent characteriz- 
ing this relaxation process was calculated [1] within 
a renormalization group approach using the method 
of e-expansion in a two-loop approximation. Later, 
the nonequilibrium critical relaxation in a short time 
regime was studied within a three-dimensional Ising 
model using methods of computer simulation [3]. 
The results confirmed theoretical predictions con- 
cerning the power character of evolution of the mag- 
netization of ferromagnetic systems, but the value 
of an exponent of 6' = 0.108(2) determined from 
these simulations significantly deviated from the the- 
oretical predictions of 9' = 0.130 (obtained by di- 
rect substitution of the parameter e = 1 for three- 
dimensional systems) and 9' = 0.138 (obtained using 
the Pade-Borel method for the summation of a very 
short series with respect to e = 1). 

According to the scaling theory, a singular part 
of the Gibbs potential $ s ; ng (t, r, /i,mo) determining 
the state of a system in the critical region is char- 



acterized by a generalized homogeneity with respect 
to the main thermodynamic variables: 



$8ing(t,T,/l,mo) 

b$ sine (b a %b a -T,b a »h,b a ™m ) 



(1) 



where t is time, r is the reduced temperature, h 
is the field, mo is the initial magnetization, b is the 
scaling factor, and aj are the scaling exponents. As a 
result, the magnetization m = —5&/5h of the system 
at the critical point (r = 0, h = 0) is characterized 
by the following time dependence: 

m(t,m ) = t-^ +1 V a <F m (m t~ a ™/ at ). (2) 

Expanding the right-hand side of Eq. @ into series 
with respect to the small parameter mot~ am ^ at , we 
obtain the following power relation: 



m(t) ~ H afc+ ° m+1)/a< ~ t e 



(3) 



All at except a m can be related to well-known crit- 
ical exponents that describe behavior of the sys- 
tem without effects related to a nonequilibrium ini- 
tial state. For this reason, Janssen et al.[2] intro- 
duced a new independent dynamic critical expo- 
nent 9'. The renormalization-group description of 
the nonequilibrium critical behavior of the system 
showed that this exponent takes positive values and, 
for t > t cr ~ to(0)~ 1 / (9 ' +,3 / 2:i/) , the initial regime 
(characterized by an increase in the magnetization 
m(t)) changes to a traditional regime of critical re- 
laxation toward the equilibrium state. The stage of 
critical relaxation is characterized by a time depen- 
dence of the magnetization according to the power 
law m ~ t~Pl vz (Fig. [TJ, where (3 and v are well- 
known static exponents determining the equilibrium 
critical behavior of the magnetization and the cor- 
relation length, and z is a dynamic exponent char- 
acterizing the critical slowing down of relaxation in 
the system. It can be shown that, when the system 
evolves from the initial ordered state with mo = 1, 
the time dependence of magnetization at the critical 
point is from the very beginning determined by the 
power law as m ~ t~^/ uz . 




FIG. 1: Schematic diagram of evolution of magnetiza- 
tion m at a critical point, starting from an initial state 
with magnetization mo. 



II. DESCRIPTION OF THE MODEL 

The critical behavior of a pure system in 
the equilibrium state is described using the 
Ginzburg-Landau-Wilson model Hamiltonian, 
which can be written as 



H GL [s] =|d d x j £ 1 [(V Sq (x)) 2 + r4(x)' 



fE4w 

\q=1 



(4) 



where s(x) is the field of the n-component order pa- 
rameter, r is the reduced temperature of the second- 
order phase transition, and g is the amplitude of in- 
teraction of order parameter fluctuations. 

Let the realization of any configuration of the or- 
der parameter in the system at a given time t be 
determined by the condition that the order param- 
eter field s(x, 0) = so(x) at the initial instant t = 
(with the initial magnetization m = mo) be char- 
acterized by the probability distribution P[so] ~ 
exp(— Hq[sq\), where 



H [s ] = / d d x — ' (s (x) - m (x)) 2 



(5) 



In the most interesting case of pure relaxation dy- 
namics of the order parameter (so-called Model A 
Pi), the exponent 9' is essentially new and cannot 
be expressed by means of the well-known static criti- 
cal exponents and the parameters of equilibrium dy- 
namics. The relaxation dynamics of the order pa- 
rameter in this case is described by the Langevin 
equation 

d tSa (x,t) = -\ SH X GL[S] +CaOM), (6) 



where Hcl[s] is the Ginzburg-Landau-Wilson 
model Hamiltonian (j4|), A is the kinetic coefficient, 
and £(#, t) is the Gaussian random-noise source, 
which describes the influence of short-lived excita- 
tions. The randomnoise source is determined by the 
following probability functional: 



P[Q ~ exp 



i I d ^ 



dt(C(x,t)f 



<Ca(M)>=0; (7) 

(Cafo t)C fi (x', t')) = 2A 8 afj 5(x - x')6(t - t'). 

Within the framework of the renormalization-group 
field theory, the critical dynamics 0, Q is described 
in terms of an auxiliary field s(x) and a generat- 
ing functional for the dynamic correlation functions 
and response functions. This functional is defined 
as follows: 

W[h,h] =lni fv{s,is)exp(-C[s,S} - H [s }) 



exp d d x dt^2(h a s a + h a s a ) \ >, 



(8) 



where C is the action functional expressed as 

s a + X(t — V 2 )s Q 



n ( 

C{s,§}= fdt fd d x J2\s a 
o a=1 [ 



0=1 



(9) 



An analysis of the Gaussian component of functional 
<j9j) for g = allows the following expressions for 
the bare response function Gq(p, t — t') and the bare 

correlation function Cq (p, t, t') to be obtained for 
the Dirichlet boundary condition (tq = oo) [2]: 

G (p,t-t') = exp[-A(p 2 +r)|t-t'|], (10) 

c ( D) ( P ,t,t') = c^(p,t-t')+c^{p,t+t'm) 

where 



C^(p,t-t') = 



1 



C^(p,t + t') 



p z + T 
1 



-A(p 2 +T)|t-t' 



p z + T 



-X( P 2 +T)(t+t') 



(12) 
• (13) 



III. RENORMALIZATION-GROUP 
ANALYSIS OF THE MODEL 

In the renormalization-group analysis of the model 
with allowance for the interaction of critical fluctua- 
tions in the order parameter, singularities appearing 
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in the dynamic correlation functions and response 
functions in the limit as r -> were eliminated us- 
ing the procedure of dimensional regularization and 
the scheme of minimum substraction [6] followed 
by reparametrization of the Hamiltonian parameters 
and by multiplicative field renormalization in func- 
tional (0) as follows: 



(14) 



where e = 4 — d, and fj, is a dimensional parame- 
ter. Calculation of all renormalization constants Zi 
(except for Zq) was described in J4|. A scheme for 
the calculation of Zq and the results of calculations 
in a two- loop approximation were presented in |2J. 
In this study, Zq was calculated in the subsequent 
three-loop approximation of the renormalization- 
group field theory. 

Introduction of the initial conditions of type ([5]) 
into the theory makes it necessary renormalize the 
response function (s(p, t)so(— p, 0)), which deter- 
mines the influence of the initial state of the system 
on its relaxation dynamics. The correction terms 
in the self-energy part of the response function ap- 
pear due to the interaction of fluctuations of the 
order parameter and are characterized by reducible 
dynamic Feynman diagrams, since they are calcu- 
lated using correlator (fTTj) . which does not possess 
the property of translational invariance with respect 
to time. Janssen et al. [2J introduced the following 
representation for this response function: 



<i) 



Gi*i(p,*) = (*(p,*)2o(-P,0)> 



= /dt'G M (p,M0lt o b,i')[ 



(15) 



In a three-loop approximation, the one-particle 
vertex function T^ (p, i)rg i with a single field in- 
sertion so is described by the diagrams presented in 
Fig. [51 which obey the requirement of containing at 
least a single Cg correlator. The factor Gi i(p, t, t') 
is determined by the equilibrium component of the 
correlator Cq in Eq. ([11]) . It should be noted that 
this correlator differs from the equilibrium response 
function G^ * (p, t — t') because the integration in 
(fT5")) with respect to time is performed starting from 
t = instead of t = — oc. However, it is possi- 
ble to establish a functional relationship [7] between 
these functions by using the functional Hql[so} (01) 
instead of ([S]) with a new interaction vertex in action 
functional ©: 



— / dt / d d X ^ (sOaSQa) (S p) 



(16) 



The additional vertex function I^q , which is lo- 
calized on the surface t = 0, appears due to averag- 
ing over the initial fields. By analogy with represen- 
tation ([T5]) . we also obtain 



-f(eg) , 



>'\ri e i)( m +"\ 



G™(p,t-t')= / dt"G 1A (p,t,t")r^'(P,t") [m] . 

Jt' 

(17) 
Solving the integral equation 

s(t-t')= f t dt"K( q ,t",t')r[^( qi t) ls{t „ )] (18) 

Jt' 

in each order of the theory, we obtain the kernel 
K(q, t" ', £'). Note that fluctuation corrections to this 
kernel in the model under consideration appear only 
in the third order (Fig. [3]). Using Eqs. $To§ and (flT[) , 
performing renormalization of the fields according to 
Eq. (I14p , and obeying the requirement of eliminating 
poles with respect to e in each order of the theory so 
that the renormalization constant Zq would remain 
finite in the limit as e — > 0, we eventually obtain the 
following expression for this constant: 



,-1/2 



dte 



'ri,o(« = o, *)[«„], 



(19) 



Upon sequential realization of the procedure de- 
scribed above and calculation of the diagrams using 
the method of e-expansion, the renormalization con- 
stant in the three-loop approximation is as follows: 



Z = l + 



n + 2 q n + 2 

-H 

6 e 12 

9\ 2 n + 2 

216 



n + 5 /, „ 1 

h e In 2 

3 I 2 



(n + 5)(n + 6) 



+ ( ^(-19 + 91n2)(n + 2) + ^(-7 + 61n2)(n + 8) Je 



1.18679(71 + 3. 13882)e 2 



(20) 



The invariance with respect to the renormaliza- 
tiongroup transformations of the generalized con- 
nected Green's function 



N,N 



([snS\«[5 ] M ) 



can be expressed in terms of the renormalization- 
group Callan-Symanzik differential equation [2J, i6j : 



N . N, 
~2 



/j9 m +CA9 a + kt8 t + I3d g + —7 + — 7 



M 



(21) 



a,/9=l 



+ Y (7 + 7o) + Cr c 7 1 a T - 1 G^ = 0. 









FIG. 2: Diagrams determining contributions to the vertex functions F\ . Curves without and with arrows correspond 
to the bare correlators and the bare response functions, respectively. The vertical line corresponds to the surface 
i = 0. 






FIG. 3: Diagrams determining contributions to the vertex functions T\ q 



(eq) 



The renormalization-group functions representing 
the coefficients in Eq. (|21[) are given by the following 
expressions: 



C= (M^)olnA = 2(7-7)) K = (^)olnr, 



7o = {^d^olnZo, 
7 = {vdfj,) \nZ s , 



7= (nd^) lnZ s , 
(22) 



where (8^)0 = ( -§- ) denotes differentiation with 

constant initial parameters g, A and r. For a 
short time regime of nonequilibrium critical re- 
laxation, the only essentially new quantity is the 
renormalization-group function 70 . In the three- loop 
approximation used in this study, this function is ex- 
pressed as follows: 



70 = -^.9 (l+(ln2-i,,, 



- 0.0988989 (n + 3. 13882) g 2 +0(g 4 ). 



(23) 



dynamic critical exponent 



e 2 ( 4 
z = 2 + - 61n--l 



(n + 8) 2 



l + e f 6 ^ + "> - 0.4384812) 



(26) 



represents a generalization of the previous results 
[9j to the case of systems with n-component order 
parameters. The final expression for the critical ex- 
ponent 9' is as follows: 



(n + 2) 
4(n + 8) 



7.2985 



el 



(it 



(n + 8) 2 



n + 3 + (n + 8) In • 




2 (n 3 + 17.3118n 2 + 153.2670n (27) 



0(e 4 



IV. ANALYSIS OF RESULTS. 
CONCLUSIONS 



The fixed point g* of the renormalization-group 
transformation is determined from the equation 
(3(g*) = 0. The general solution of differential 
Eq. (f2"Tj) by the method of characteristics at the fixed 
point has the following scaling form [2] 



G^({a;,i},r,T - 1 ,A, 5 *, Ai ) = ^- 2+ ^f 
x G^({lx,l 2+ <*t},rr 2 +** ,r^l 2 +<' ,\,g*,n), 



(24) 



The series of e-expansions exhibit factorial diver- 
gence, but they can be considered in an asymptotic 
context [11| • In order to obtain physically reasonable 
values of the critical exponents for three-dimensional 
systems at e = 1, special methods for the summation 
of asymptotic series have been developed (12h17|. 
the most effective being the Pade-Borel, Pade-Borel- 
Leroy, and conformal mapping techniques. We used 
the Pade-Borel method to perform the asymptotic 
summation of e-expansion series (|27p for the critical 
exponent 8'. According to this method, the series 



where rj s — 7*, rjs = 7* and 770 = 7o ~ are the 
exponents of anomalous dimensions. The functions 
entering into Eq. (|24|) can be related to the criti- 
cal exponents involved in the scaling relations, for 
example: 



00 

n=l 



is replaced by the integral 



(28) 



z = 2 + C, l/v = 2-n*, 



7o 



(25) 



C*+7* + y )/(2 + C), 



which determine the critical relaxation dynamics 
(z), correlation length (is), and nonequilibrium evo- 
lution (6') of the magnetization. In determining the 
h value, we used the data of Kleinert et al. [8( on 
the coordinate of the stable fixed point g* and the 
results of calculations [9] of the dynamic critical ex- 
ponent z for the Ising model, which refined (in the 
three-loop approximation) the previous values of ex- 
ponents for this model. The final expression for the 



'(e) = / e- t B{et)dt, 






(29) 



where B(x) is the so-called Borel image. In con- 
trast to the initial series (|2"5)) having a zero radius 
of convergence, the Borel image has a finite radius 
of convergence determined by the parameters of the 
asymptotic nth term of series (|29|) for n ^ 1. Then, 
the Borel image is subjected to the Pade approxi- 
mation, according to which B{x) is replaced by a 



TABLE I: Results of calculations of the critical exponent 
9' for the Ising model and XY model in comparison to 
the results of computer simulations 



Calculation method 



Critical exponent 9 
Ising model XY-model 



Two-loop approximation 



Substitution s = 1 
Pade-Borel summation 



0.130 
0.138 



0.154 
0.170 



Three-loop approximation 



Substitution s = 1 
Pade-Borel summation 



Computer simulation 



0.0791 
0.1078(22) 



0.0983 
0.1289(23) 



0.108(2) [3] 0.144(10)[18] 



rational function of the following type: 



L 
J2 a i x% 

i=0 



[L/M] = M 



(M > 1), 



(30) 



j=o 



the expansion of which into Taylor's series (in the 
vicinity of X = 0) coincides with that of the Borel 
image as far as possible. The function of type (|30|) 
has L + 1 coefficients in the numerator and M + 1 
coefficients in the denominator. The entire set of 
coefficients is determined to within a constant fac- 
tor (for certainty, bo = 1), so that there are a total 
of L + M + 1 free parameters. This implies that, 
in the general case, the coefficients of expansion of 
the [L/M] function into Taylor's series must coincide 
with the corresponding coefficients of series (T2"5]) . If 
this series has a finite number N of terms, the num- 
ber of coefficients in expansion of the [L/M] function 
must obey the condition L + M < N. 

In the case of a three-term series under considera- 
tion, the summation of the e-expansion for exponent 
9' was performed in the [2/1] approximation. In the 
given order of the theory, the choice of this approxi- 



mation is (in accordance with the results of analysis 
performed in [17]) preferred for obtaining more ac- 
curate value of the sum. The results of calculations 
are presented in the Table [TJ 

A comparison of the results of our calculations of 
the critical exponent h to the values obtained by nu- 
merical modeling within the Ising model [3] and XY 
model [18] using the method of short time dynamics 
(see Table|l| clearly demonstrates that the values ob- 
tained in the three-loop approximation better with 
the results of computer simulations as compared to 
the results of a two-loop approximation. 

In this study, we have presented a field theory 
description of the nonequilibrium critical relaxation 
of a system within the most interesting dynami- 
cal model A (according to the Hohenberg-Halperin 
classification [l|). It is shown that only beginning 
with a three-loop approximation does the theory of 
these processes involve an additional vertex func- 



tion r 



(eq) 

1,0 



localized on the surface of initial states 



(t = 0), which provides fluctuation corrections to 
the dynamic response function due to the influence 
of nonequilibrium initial states. As a result, only 
allowance for these fluctuation corrections (reflect- 
ing the influence of the nonequilibrium initial states) 
ensures adequate description of the relaxation pro- 
cess. Using this three-loop approximation and the 
method of e-expansion, it is possible to obtain the 
values of independent dynamic critical exponent 6' 
describing the evolution of the system during a short 
time evolution in close agreement with the results of 
computer simulations. 
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